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I, INTRODUCTION 


Control system complexity and demands on reliability are 
rising with ever increasing automation. More and more control system 
engineers are looking to solid-state digital techniques for new 
answers to the old problems of _ reliability, complexity, and 
economy, Digital techniques provide many advantages such as accuracy 
and computational speed as well as flexibility and versatility, 
Which analog techniques could not provide, By the use of digital 
techniques in control systems, sensitive, and efficient control 
devices can easily be built, Another’ significant advantace of 
digital control systems is that the utilizing of digital components 
in a control system allows time sharing of important parts 
of the system which resuits in economy in the use of equipment. 
The digital elements also require much less power than the analog 
elements, The data signals in digital control systems can also be 
transmitted and received in pulse-code form almost’_ error-free, 
even in transmission through a noisy media, Furthermore, 4 
digital computer in control gives another important adventage, it 
makes feasible system compensation by nonlinear programming and 
adaptive or self-optimizinzg control. 

The digital servo system is becoming widely accepted in 
machine tool control, radar systems, navigationel system, and 
many other fields. A digital servomechanism uses numerical 
information to provide a simple but highly accurate method of 
positioning control, Figure 1 shows a block diagram of a typical 


digital servo system. The major components of the basic digital 
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servo system are: 

1, Controlled system 

2, Control elements 

3. Analog-to-Digital converter 

4, Digital-to-Analog converter 

5. Digital data processing equipment or computer 
The controlled system is normally a linear servomotor but step 
motors have been used for quite some time in both closed loop and 
open loop digital positioning systems, Tne digital computer is used 
as the eenbres control unit for data processing, system compensa- 
tion, and decision making in some cases, In a Simole digital servo 
System, a Berea racer such as an up-down counter can be used. 
The digital. toeanalog and analogeto-digital converters are used to 
interconnect the digital and analog quantities fon controlied and 
control elements. One familier type of analog-to-digital converter. 
encoder is the shaft-position encoder disk employing a gray or 
binary code, This digitizer is used to translate .the shaft 
position into digital notation by making use of a _ mechanically 
rotating shaft with coding disk centered on the shaft. Figure 2 
illustrates a typical coding disk for a shaft digitizer, The 


brushes contact the disk along the reading line at the point 
shown with a single brush positioned “in : the center ‘of each 


- Be ener, 


concentric band, These concentric rings represent a 2 
and ‘each . ring.: is °. constructed of several segments made of 
either conducting material (the darkened areas) or “insulating 


material ( the wunshaded areas), A signal is connected to the 





conductor, and if a given brush makes contact with a segment of 
conducting material, a "i" signal will result, but if the brush 
is over the insulating area, the output from the brush will be 
a "0", The resolution of this encoder is determined by the 
number of channels, For a 4echannel coding disk, a resolution of 
one part in sixteen is achieved, Two types of encoder disk, the 
gray and binary code disks are shown in Figure 3, The coded disk 
may be made of transparent and opaque materials arranged in coded 
patterns, The photocells and sources are used to read the angular 
position in a digital code, 

In operation, the analog-toedigital converter converts the 
analoz output of the controlled system into a-numerical quantity, 
which is compared with a numerical order, and the difference is 
converted by the digital-to-analozg converter into an analog signal 


capable of driving a servomotor until the difference is eliminated, . 


A, OBJECTIVES 

The methods for analysis and synthesis of the analog system 
such as Bode plots, Root Locus, Parameter plane plots, etc, have 
been established for a long time, These methods work very well, 
however, for the discrete or digital system, there is no such 
method that is as good as those for the analog = system, Even 
the well known z-transform method can become very complicated if 
the order of the systems analyzed are higher than third or fourth 
order, 


The objective of this paper ‘was to investigate a typical 


16 





digital positioning servo system in order to obtain some good 
methods for analysis and synthesis, In the proposed system, light 
travels from a light source through equally spaced openings in the 
code disk to activate the sensor. The sensor, in turn, generates 
feedback pulses at equilibrium angular position. The sampling 
period of this system is not constant, depending on the velocity 
of the motor, the z-transform method of analysis and design cannot 
be used here, Another objective was to find a suitable mathematic 
representation of such a_ system, 

All of the studies had been done by simulation using Digital 
Simulation language DSL/360 on the IBM 360 computer system, this 
language made the simulation quite easy, however, care must be 
taken in order to obtain accurate integration in the system with 
high velocity. 

In Chapter II, the system identification and method of 
Simulation are given, Chapter III presents the discussion of the 
results obtained from the simulation of second and third order 
servo systems, This includes the study of the stability limit, 
percent overshoot, and the effects= of “variation .-of the 
servoamplifier gain. This chapter also shows the results obtained 
from investigating the behavior of the system with ramp input, 
The Bode plot and parameter plane plot were used to compensate 
the analog system in Chapter IV. The results were then used in 
simulating the digital system, since there was no direct ‘compensation 
method, 


The results were obtained both in graphs and printed output, 
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however, only graphs are included, Listing of all computer programs 
are presented after the conclusion and recommendation in Chapter V, 

In each study, the digital system with several numbers of 
opening slots and the continuous counterpart were simulated, The 
studies of the digital systems were based on the results of studies 


of the continuous system, 
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Il, TH DIGITAL POSITIONING SERVO 
A, PHYSICAL PROPERTISS 

The digital positioning servo system under consideration consists 
ou ame digital comparator, a digital-toeanalog converter, a controlled 
system, and an analogetoedigital converter, The analog-toedigital 
converter consists of a low inertia encoder disk mounted tightly 
on the servomotor shaft, two photoelectric sensor units, and a 
direction-sensing logic circuit, The encoder disk has uniformly 
spaced slots around its circumference, Each photoelectric sensor 
unit consists of a light source and a photodiode, Figure 4 
illustrates the encoder disk and the sensor units and Figure 5 
illustrates the direction-sensing circuit. 

Each time an opening passes through the beam from the light 
Source, 4 pulse is generated by the sensing unit, and it will 
eventually be fed to the upedown counter which acts = as the 
digital comparator, The second sensor unit is for the purpose of 
direction sensing, When the servomotor turns in forward direction, 
the feedback pulse becomes a downeenable signal, when the servomotor 
turns backward, the feedback pulse becomes an up-enable signal. 
The two sensors are adjusted so that their respective outputs, X* 
and Y, overlap as shown in Figure 6, Figure 6-A and 6-B show 
the inputs and outputs of the direction-sensing circuit when the 
servomotor turns forward and backward, respectively. Figure 7 shows 
a typical digital-to-analog converter, 

For an input as a step function, an appropriate binary number 


representing the size of the step function is fed into the up-down 





counter, The output of the up-down counter is, in turn, cenverted 
into an analog signal and is sent to the  servoamplifier to 
drive the servomotor, The sensing device will measure the ‘shaft 
position, encode, and send this information to the up-down counter, 
The positioning process will be finished when the binary number 


in the up-down counter becomes zero, 


B. COMPUTER SIMULATION 

The DSL/350 language was specially designed PoE use in 
simulation. DSL/360 provides many advantages over other languages 
and subroutines in existence, especially in simplicity and time 
saving. DSL/360 not only provides these advantages in simulating 
the linear physical system, but in simulation of the non-linear 
system also, L/360 provides a basic set of function blocks 
from which a physical system may be modeled such as integrators, 
quantizer, pulse generators, function generators, and so_ on, In 
addition, FORTRAN Library functions from the Scientific Subroutine 
Package may be utilized, The user may also provide his own 
subroutines if necessary. 

In programming, first, each part of the digital positioning 
servo system was separated into single function biocks as shown 
in Figure 8, Next, each single function block was replaced by a 
related DSL function or standard FORTRAN Library function as 
appropriate, Finally, all the DSL and FORTRAN Library function were 
put together:'in the DERIVATIVE REGION of the program, One advantage 


obtained in using the DERIVATIVE REGION is that the functions are 





not required to be in any specific order, 
In the early parts of simulation, The servomotor was assumed 


to be 2 second order motor with the Transfer function 


K 
v 


Gils) = ——— . 
a s(s +1) 


Later, the third order servomotor with the transfer function 


KY 


Cee EE 
as s(s + 1)(0.2s + 1) 


was used, The FORTRAN Library function 


Y = AINT(X) 
was used to modélled the quantizer. This function was also used 
to quantize the ramp input function. The gain of the digital. 
tomwanalog converter was designed to be equal to 1/K, to make the 
system a unity feedback system. 

The fourth order Runge-Kutta with fixed integration interval 
(INTHG RKSFX) Integration routine was used in all simulations to 
assure the accuracy of the results, It is very important to 
choose a correct integration interval (DELT), since too large DELT 
will cause inaccurate integration and too sm ll DELT will result 
in excessive computing time, The DELT used in each simulation 
is approximately five to ten percent of the time interval between 
two consecutive pulses generated by the digitizer at its maximum 
angular velocity in that specific simulation. 

In miltiple runs, the TERMINAL region was used to terminate 


the previous run, to provide the new sets of parameter values, 


2) 





and to terminate the program after the output data.of the last 
run had been plotted. The FORTRAN subroutine DRWG was used in 
the SAMPLE region to collect points for each sample interval 
(DELS) to be plotted by the CALCOMP XY Plotter, 

In most cases, The digital system with 10, 45, 90, and 180 
openings and the continuous system were simulated and the results 
were compared, A system with one opening was used in stability 
limit studies since it was the worst case, To compare the outputs 
of the digital system with output of the continuous system, the 
inputs of all systems must be equal. If a step input to the 
continuous system is 1.25664, inputs to the digital system with 
10, 45, 90, and 180 openings must be 2, 9,.18, and 36 respective- 
ly. 

Thess input values can be readily calculated by the method 
shown in example below, 

EXAMPLE 

To calculate comparable step inputs to the continuous and 
digital system with 10 openings. 

From Figure 8, 


vi =Re-C 


in(continuous ) 
R - iGBIT 
V i pe 5 eee 
in(digital) K 
2 (R, ~ IGBIT) 
OP 


for zero initial conditions, 
V , = 
in(continuous ) 


2a 





Ze 
V pee te 
in(digital) OP 
to have equal inputs to servomotors, requires 
Zk 
n 
OP 





by arbitarily choosing RO = 2 
R = 1,25664 
The values of inputs to digital system with other numbers of 


opening were also based on this’ selection, 
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it. DlsCuss iON 


The first part of this chapter contains a discussion of the 
response of the second order servo system to a_ step function 
input, The third order system with both step and ramp function 


input is discussed in the later part of the chapter, 


A, STEP RESPONSE OF THE SHCOND ORDER SYSTEM 

The main concern in the response of the second order servo 
system was its stability limit. The forward gain of the system 
had been increased until the stability limit was obtained, 

The model of the second order continuous system used in this 
Simulation is always stable regardless of the size of the gain, 
The digitel system was also stable at all value of gain regardless 
of the number of openings in the position encoder, The gain 
of the system was increased to 10° to make sure that the system 
would always be stable, Figure 9, 10, 11, 12, 13, and 14 show the 
step response and plots of the binary numbers in the up-down 
counter versus time for the digital system with number of openings 
equal to i, i0, 45, 90, and 180, and for the continuous system, res- 
pectively. Figure 15, 16, 17, 18, 19, and 20 show their respective 
phase plane plots. The plots of the output of the up-down counter 
were not scaled correctly, however, they could be read easily. 
Each quantizied step of the Ne plot represents a one, 

The difference that can readily be seen in the responses of 
the digital and continuous systems is the magnitude of both the 


angular displacement and velocity. The difference is very pronounced 
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when the number of openings is small, The digital system with 
fewer openings has longer time delay between two consecutive 
feedback pulses, This time delay has allowed the motor to move 
farther with higher velocity, similar to the response of the 
system without a feedback loop, before the excitation signal is 
reduced by the next feedback pulse, The system with a large 
number oof openings has a very high sampling rate, therefore this 
system behavior is much closer to that of the continuous. system, 
a code disk with 20000 openings has been developed. This disk 
gives the system almost the same response as that of a continuous 
system, 

The phase plane plots of the digital system are not as 
smooth as those of the continuous system. They actually consist 
of many short sections, tach portion represents the delay 
interv2l between two consacutive pulses, This can easily be observed 
in the phase plane plot of the system with lower number of 
openings such as shown in Figure 16, 

Figure 21. contains the comparison of the angular displacement 
of the analog and digital system with 10 and 45 openings. This 
plot clear]: shows that the natural frequency of this second 
order digital system is smller than that of the analog system, 
It indicates that the root locus of the digital system would 
cross the W-axis closer to the origin than the root locus of 
the continuous system, This distance becomes smaller if the number 


of openings is decreased, 
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The plots of the step response of the digital system also 
consist of many short sections, It is very easy to observe this 
phenomenon in the response of the system with very low. gain 
such as those shown in Figure 22 and 23, Figure 24 demonstrates 
how this phenomenon occurs, At the beginning of the simulation, 
the system behaves like a system without a feedback loop. The 
rate of increase of velocity is very high at first, however, it 
becomes lower with time, The displacement plot also becomes more 
linear with time, When the displacement is large enough, a feedback 
pulse is generated, the output of the up-down counter is changed 
by one, the velocity and the rate of change of the displacement 
are changed to follow the new input voltage to the servoamplifier. 
This process occurs repeatly until the output of the up-down 
counter is eliminated, and the system eventually comes to rest, 

Figure 22 and 23 also illustrate that the digital system behaves 
like the continuous system in another aspect. When the gain is 
reduced, the peak overshoct decreases, This peak overshoot can be 
eliminated by making the gain small enough. It is also observed 
that, for all cases studied, the steady state displacement of the 
digital system is always larger than that of the continuous 
system with comparable size of step input. The difference, however, 
has never exceeded the angular displacement between two adjacent 
openings, The system with more openings in the encoder disk gives 
smaller error than the system with fewer openings. The maximum 
error for the system with 10, 45, 90, and 180 openings are 


less than 0.6283, 0.1396, 0.0698, and 0,0349 radians respectively. 
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B, THIRD ORDER SYSTSM 
A third order system was studied in greater detail than 
the second order’ system, In addition to what had been studied 
from the second order system characteristics, the variation of 
the percent overshoct was investigated in this part, 
1, Stability limit 

A third order continuous system with the transfer function 
Shown in Chapter II reaches the stability limit when the 
Servoamplifier gain, K, is equal to 6, The digital ~ systems, 
however, require larger gains to reach their stability limits, 
Figure 25, 26, 27, 28, and 29 illustrate the phase plane plots of the 
third order digital system with the gain equal to 6 and the 
number of openings equal to 1, 10, 45, 90, and 180 respectively, 
Figure 30 shows the phase plane plot of the continuous’ system 
with the same gain. The phase plane plots of the third order 
digital system converge into a limit cycle and can clearly be 
observed in the case the system has i or 10 openings. By 
inspecting the characteristic of the phase plane plots, especially 
in the case where the number of openings is small, it is 
obvious that the third order digital system requires higher gain 
to reach the stability limit than the continuous system, The 
time delay also causes the displacement, velocity, and limit 
cycle size of the digital system to be larger than those of tho 
continuous system, 

Since there had not been any known method that can be 


used to find the gain at stability limit of this digital servo 
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‘System, a trial and error method was used to find this gain, the 
gain of the third order digital system with 10 openings at its 
stability limit was found to be approximately 6.25. This gain 
would be reduced toward 6 if the number of openings is increased, 
The size of the limit cycle is also decreased with the increase 
of the number of openings. 

The results from simulating the third order digital servo 
system with number of opening equal to 10 as illustrated in 
Figure 31, 26, 32, and 33 also show that the size of the limit cycle 
increases with the gain. The increase of the gain also causes 
the system to reach this limit cycle faster, Upon making the 
gain greater than the stability limit gain, The phase plane plot 
diverges for a period of time before reaching another limit cycle 
as shown in Figure 34, For higher gain, The phase plane plot 
diverges into a limit cycle and stays in the limit cycle for a 
period of time then diverges to infinity as shown in Figure 35, 
With the gain a little higher as shown in Figure 36, the phase 
plane plot diverges to infinity but has a tendency to go into 
a limit cycle at first. When the gain is much higher, the 
system just becomes unstable as shown in Figure 37, 

2. Percent overshoot 

For a continuous and stable servo system, the overshoots 
of the ste n response always decrease with time, and the plot of 
the log of the percent overshoot versus time is a straight 
line, The step response and the plot of log of percent overshoot 


of the continuous system are shown in Figure 38 and 43 respectively, 
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To investigate the behavior of the overshoot of the digital 
positoning servo system, a gain of 2.5 was used. With this gain, 
¢ and Wn of the continuous system are equal to 0.166 and 1,54 
respectively, Figure 39, 4, 41, and 42 illustrate the step response 
and output of the up-down counter of the third order digital servo 
system with 180, 90, 45, and 10 openings respectively, The step response 
of the continuous system shows 6 prominent overshoots, the responses 
of the digital system with 180, 90, 45, and 10 openings show 5,,4, 
3, and 2 overshoots respectively, The time at which each of the 
corresponding peaks of all systems occur is the same, Figure 43 
illustrates the plots of the percent overshoot of the responses 

f the continuous and digital system with 180, 90, 45, and 10 
openings versus time, The numerical values of percent overshoots 
are included in table I. The inputs to all systems were adjusted 
so that they would be equal and the percent overshoot could be 
calculated from the same basis, the steady state output of the 
continuous systen, 

The overshoots of all digital systems are higher then the 
corresponding peaks of the continuous’ system, The difference, 
however, does not exceed the angle of the two adjacent openings, 
The digital system with smaller number of openings has_ larger 
percent overshoot than the system with larger number of openings. 
The steady state response of the digital system is also larger 
than that of the continuous system. The difference is also less 
than the angular displacement between two adjacent openings. 


At the beginning, the plot of the percent overshoot of the 
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digital system tends to be parallel to that of the continuous 
system. The plot of the percent overshoot then diverges to the 
smallest peak which is slightly larger than the angular displacement 
between two adjacent openings, Since’ the system with fewer 
openings has larger displacement between the openings, its response 
requires less time to reach the last overshcot,. The number of its 
overshoots is also less than that of the system with more openings, 

The responses of all the digital systems that have less 
than 10 openings were also studied, All of the responses heve 
two overshoots and the characteristics are similar to those discussed 
above, However, only the response of the system with one opening, 
which is the extreme case, is illustrated in Figure 44, 

3. Resvonse to ram input 

The ramp input function, Bt, was also applied to the third 
crder digital positioning servo system in order to investigate its 
stability. The servoamplifier gain used in this study was 2.5 
and the slope of the ramp was varied from 0,5 .to .5, . Practically, 
a method of applying a ramp input to the digital servo system 
is to apply a series of impulses into the up-down counter, The 
slope of the ramp can be controlled by varying the period of the 
impulse, A digitized ramp function was used in this simulation 
instead of the impulses. The study of the response of the digital 
servo system with ramp input was done by simulating the system 
with 10 openings only. It had been shown from the previous studies 
that the digital system with more openings behaves closer to the 


continuous system than the system with fower openings, 
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Figure 45 and 46 illustrate the response of the digital system 
with 10 openings and the slope of the ramp input equal to 0,5, 
PO, 4a, 250, 2.5, 3.0, 4.0, and 5,0. When the slope is small, 
for example B= 0.5, the rate of the input pulse applied to the 
up-down counter is slow and causes the velocity of the shaft and 
the feedback sampling rate to be too slow and nonuniform, The 
nonuniform and slow rate of output of the up-down counter, in 
turn, cause the shaft to turn abruptly and the system becomes 
unstable, The longer time delay is also an important cause for 
‘the instability of the system, The plot of velocity and output 
of the up-down counter for the digital system with the slope of 
the ramp input to the system equal to 0.5, 1.0, 1.5, 2.0, 2.5, 
3.0, 4.0, and 5,0 are shown in Figure 47, 48, 49, 50, pears - 
and 54 respectively. when the slope of the ramp input is increased, 
the input pulses are fed to the up-down counter at higher rate, 
the servo shaft turns faster, the feedback sampling rate beccmes 
higher, the output of the up-down counter becomes more linear and 
the velocity of the shaft becomes more stable, Unlike the continucus 
system, the velocity of the digital system will not be constant 
no matter how large a slope is used for the ramp input, It rather 
varies within a small limit, This limit cycle is caused by the 
time delay between pulses, The limit cycle will be smaller if 
the time dolay is. shorter, 

Another interesting result observed ‘is that the rate of 
decrease of the magnitude of the limit cycle with the increase 


of the slope of the ramp is nonuniform, At first the magnitude 


ap 





of the limit cycle becomes very small when the slope is increased 
to 2.5. However, when the slope is increased to 3 and 4, the 
magnitudes of the limit cycles are some what larger, The limit 
cycle is reduced to a smaller size again when the  silope is 


increased to 5, 
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IV. SYSTEM COMPENSATION 


Since there had not been any known method that can be used 
to compensate this digital positioning servo system directly, the 
author chose to compensate the comparable continuous system with 
known methods and simulated the compensated. system as both a 
continuous and a digital system, The results from the simulation 
were then compared and studied, The continuous and the digital 
system with 10, 45, and 180 openings were similated and their 
step responses were plotted on the same graph to make the 
comparison more obvious, The simulation of the digital system with 
90 openings was omitted to prevent the graph from being overcrowded, 
The system compensated was the third order servo with servoamplifier 


gain equal to 6, 


A, LAG COMPENSATION 

Since the programs and computer facility were available, the 
parameter plane plot method was used in lag-compensating of the 
system. This method is very efficient and flexible and the computer 
make it very easy to obtain the plot, The author selected the 
PARAM M program to use in this compensation, The complete details 
of this program are shown in the appendix. 

In order to obtain the parameter plane curve, the characteristic 
equition of the system and other requirements such as the value 
of ie, Gin, , or OG needed to be plotted and the scales of the 
plot must be provided as the input data to the PARAM M program, 


Figure 55 illustrates the block diagram of the third order 
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continuous system with the lag compensation section. The lag 
compensation block has the transfer function 
(Ts + 4) 


G (s) = 
(Ts + 1) 


The transfer function of the third order system with the lag 
compensation becomes 
Ts + 1) 6 


( 
G(s) =—_—. 
- (Ts 1) sis + 1)(@ees 4.49 





or 


_p (s +2) 30 
2(s +p) s(s+1)(s + 5) 


The characteristic equation of this system is 
s+ (6 + ps? + (5 + 6p)s* + (5p + 305 )s + 30p ='0 


let c= 


N | ‘OC ae) 


6 = 


the characteristic equation becomes 


s’ + (6 4a)s> + (54+ 6%)5° + (5 +30B)5 + 30a= 0 


The parameter plane plot of this characteristic equation is shown 
a7 Figure 56. am and @B are the unknown pair of variables which 
their values can be directly selected from the suitable value of 
S$ and Wy. The values of pole and zero can also be readily 
obtained from the values of a@ and B. Three pairs of poles and 
zeroes were chosen and simulated, The step responses of the 


compensated continuous and digital system with 10, 45, and 180 





openings are illustrated in Figure 57, 58, and 59. The values of 


pole, zero, €, @n, and predicted and resulted values of Mot i 
and te of the compensated system are included in table II, 

The simulated results show that the characteristics of the 
digital system are similar to those of the: continuous system, 
The results also agree with those studied in the previous chapter, 
The magnitude of the peak overshoot and steady state response of 
the digital system are always larger than those of continuous 
system, All the differences are also within the limit of the 
angular displacement between two adjacent openings, 

Care must be taken in selecting the value of the zero for 
the lag compensation section, If this zero was too close to the 
origin, the real root would dominate the system instead of the 
selected pair: of imaginary roots, If the zero was too fer 


from the origin, the value of § wovld be too smaji and cause 


the system to be very slow, 


B, LEAD COMPENSATION 

The same characteristic equation was used to obtain tne parameter 
plane plot for lead compensation, However, the scales, the starting value 
of Wy, the required values of §& and Wn curves, and etc, must 
be reselected in order to move the curves into required region 
for the lead compensation, The block diagram of the system with 
lead compensation section is the same as that of the lag 
compensated system, Its parameter plane plot is shown in Figure 60. 


The pole and zero of the lead compensation section were kept 
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within one decade apart to prevent the severe reduction of the 
forward transfer gain of the section, This imposes another limitation 
into the selection of the value of 6. The value of £@ _ should 
not be larger than 10, Three pairs of poles and zeroes’ were 
selected and simulated, Figure 61, 62, and 63 illustrate the step 
response of the compensated system in both continuous and digital 
versions, Table III consists of the values of pole, zero, 0, Gn 


predicted and resulted values of “ae ee and t. of the compensated 


system, Since the natural frequency of the lead  compensatéd 
system is larger than that of the lag compensated systen it 
results in smaller and t.. The other results of this = study 


are Similar to those described in the lag compensation study above. 


C. LAG-LEAD COMPENSATION 

The block diagram of the third order continuous system with 
lag and lead compensation sections is shown in Figure 64, There 
are four variables in this system which make it more difficult 
to use the parameter plane plot method in compensation, The Bode 
plot method was used in this compensation. This method is rather 
simple but requires some degree of neatness. Figure 65 and 66 
illustrate the Bode plots of two different lag-lead compensations. 
The limitation in selecting the values of the zero of the lag 
section and the spread between the pole and the zero of the lead 
section are also imposed here, Both of the compensated systems 
were also simulated and studied, The results of the simulation 


of the lag-lead compensated systems are shown in Figure 67 and 68. 
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Table IV consists of the values of pole, zero, f, Wn; predicted 
and resulted value of ey ote and te. The results of the lag. 
lead compensation study are also similar to those of lag and lead 
compensation studies, 
D, COMPENSATION DISCUSSION 

The results from all compensation studies as well as those from 
the previous chapter show that the steady state response of the 
digital system is always within the angular displacement between 
two adjacent openings larger than that of the comparable analog 
system, However, this difference is not constant, it changes randomly 
within the limit, The exact position of the steady state response 
of the system can not be determined. It seems that the best method 
to compensate the digital system is to compensate its analog 
counterpart, 

The maximum error of the system under investication is found 
to be 


1.296 x 10° 


OP 


seconds 


In comparison,. the analog system has the accuracy of approximately 
2 minutes of arc and the digital system with gray or binary 


shaft-position encoder has the accuracy of 


6 
1,296 X 10 seconds 

an 
where n is the number of the conducting channel (see ref, 5). In 
order to have the same accuracy as of the analog system, this 


digital system must have approximately 10000 openings in the encoder 
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dusk, The digital system with gray or binary code disk requires 
13 channels in the code disk to acquire comparable accuracy, It 
is trivial that this digital system has an economic advantage over 
the gray or binary code disk system, This system requires only 
two sensing units for any number of openings, The gray or binary 
Sshafteposition enceder unit requires one sensing device for each 
channel in the code disk and the sensing circuit has to be more 
complicated, 

Another adventage is that all the studies show that the 
response of this digital system is always larger than the response 
of the continuous system, This observation could be used to reduced 
the ~steady state error by making the size of the step input to 
be less than the required oe to the system, The difference 
should be as close to a half of the maximum error as _ possible, 
The difference of exactly a half of the = anguiar displacement 
between two adjacent openings would reduced the steady state error 


to within one half of the maximum error, 
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CONCLUSS ION 


The characteristics of the digital positioning system under 
investigation can be very similer to the analog system if the 
feedback sampling rate is high enough, This system has some 
characteristics which are superior to the conventional'analog and 
Sampled-data system, especially in accuracy and _ economy, 

Compensation of this system can be accomplished by using the 
results obtained from compensating its analog counterpart by 
conventional methods, The accuracy can be readily improved by 
modifying the encoder disk. 

From the discussion above, it can be seen that there are 
many areas left to be studied, The author believes’ that 
further study of the mathematical analysis and of the direct method 


of compensation are necessary, 
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Figure 1, Block diagram for a typical digital control system. 
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Figure 2, A typical shaft position-to-digital converter 
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Figure 4, Code disk and photoelectric sensor units, 


45 








*YPNOIFO SUT SUSS=“UCTYOeATG °¢ eMnsty 


TWNOIS GIG Wa-NMOG = 
© 


Ind1no ete <= LNANI 
TWNOLS FI1awa-dn a= 


al 





ATNoAFO BFuUPSsues~uoTyzoodttp eyzy Jo sTeusts qandyno pue ynduf 


paemyoeq suUmMy 10Z0U = (q PaBMIOJ SUMZ TOQOU (ze 








°g OmMsty 


"Oo 
TWNOIS 
4 STGWNE-NMOC 
. TyNOrs 
a, A TaWa-dn 
if 


LAdINO 


LAGNI 


= 4 








LNdINO 


wel 


JOYAGAUOD SOTEUS“04-TeqTZTG Veemanrs 


eer 


De eS ge oe 


waINNnOD a/n 


WOU 


y SINdNI u 


48 








a) continuous system 


a (ee cS lia = a a ee a 
v/D 


counter converter 


- 
| 
= 
r 





b) digital system 


Figure 8, Positioning serva system in 


single function block diagram 
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Figure 15, Second order digital servo, K =107, Rais OP=1 
phase plane plot 
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Figure 16, Second order digital servo, K =10", Rie, OP=10 
phase plane plot 
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Figure 17, Second order digital servo, K =10°, R=9 OP=45 
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phase plane plot 
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Figure 18, Second order digital servo, K,=10', R=18, OP=90 
phase plane plot 
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phase plane plot 
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Figure 20, Second order analog servo, K =10", R=1, 25664 
phase plane plot 
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--Figure 22, Step response of second order digital servo 
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Figure 25, Third order digital servo, K =6, R =i, OP=1 
phase plane plot 
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Figure 27, Third order digital servo, K =6, R=9, Oop=45 
phase plane plot 
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Figure 28, Third order digital servo, kK =6, R =18, OP=90 
phase plane plot 
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Figure 29, Third order digital servo, K =6, R36, OP=1 80 
phase plane plot 
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Figure 30, Third order analog servo, K =6, R=1, 25664 
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phase plane plot 
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Figure 32, Third order digital servo, K =6.1, R=2, OP=10 
phase plane plot 
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Figure 33, Third order digital servo, K =6.25, Ru=2, OP=10 
phase plane plot 
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Third order digital servo, K =6.5, R=2, OP=10 
phase plane plot 
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Figure 36, Third order digital servo, K =7.2, R=, OPp=10 
phese plane plot 
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Figure 37, Third order digital servo, K =7.5, R =e, OP=10 
phase plane plot 


78 





radians 


Displacement, 


<.0e 


87 


Le) 

nn 

fas) 

s> 

m 

oO 

I 

ed 

o-4 

ari 

I 

wt 

Ls] 

i] | 
6.0d 12.¢¢ 18.00 24.00 39 


Time, seconds 


Figure 38, Third order analog servo, K =2.5, R=1,25664 


displacement versus time 
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Figure 40, Third order digital servo, K =2.5, R=18, OP=90 


displacement and up-down counter output vs tinie 
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Figure 41, Third order digital servo, K =2.5, R=9, OP=45 


displacenent and up-down counter output vs time 
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Figure 42, Third order digital servo, K =2.5, R=, CP=10 


displacement and up-down counter output vs time 
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Figure 43, Third order servo systen, K=2.5 
OP=10, 45, 90, i180, and analog 


percent overshoot versus time 
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Figure 44, Third order digital servo, K =2,5, Ris2, Op=1 


displacement versus time 
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Yamin input, 8=0.5, 1.0, 2.5,and 2.0, Gr=10. K =2.5 


86 





radians 


Displacement, 


a 
~ 


th 


ie 


Pai OG 


Lh. OG 


oS 


re) 





ae 208 17.93 18.09 2 
Time, seconds 


Figure 46, Third order digital servo, displacemsnt vs time 
ramp input, B=2.5, 3.0, 4.0,and 5,0, ,OP=10, K =2.5 
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Figure 51, Third order digital servo, ramp input, .B=2.5, K 
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Figure 54, Third order digital servo, ramp input, B=5.0, K =2.5, OPp=10 


velocity and up-down counter output versus time 
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Figure 57, Lag compensated third order positioning servo 


p=0,013, 220,18, OP=10, 45, 180, and analog 
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Figure 58, Lag compensated third order positioning servo 


p=0.017, 2=0.19, OP=10, 45, 180, and analog 
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Figure 59, Lag compensated third order positioning servo 
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Figure 60, Parameter plane plot for lead compensation 


of the third order positioning servo 
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Figure 61, Lead compensated third order positioning servo 


p=i4.5, 2=1.85, OP=10, 45, 180, and analog 
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Figure 62, Lead compensated third order positioning servo 


Pai5,0, 2217500 CF =NGe 4S, 18@. and analog 
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Figure 63, Lead compensated third order positioning servo 


p=19,0, 2=1.90, OP=10, 45, 180, and analog 
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Figure 67, Lag-lead compensated third order servo 
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Figure 68, Lag-lead compensated third order servo 


=),2, z~. -=0.5, Py pag nies Ds er ee 


Prag lag, 
OP = 10, 45, 180, and continuous 
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APPENDIX 


This appendix presents four examples of the programs used to 
Simulate the system under investigation and a program used to plot 


the parameter plane curves, 
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TITLE SHCOND ORDER ANALOG POSITIONING SERVO 
3 STEP RESPONSS AND PHASE PLANS PLOT 


* 
* 
* 


INTHG RKSFX 
INTGER NPLOT, NCUR 
CONST NPLOT=2, NCUR=2,KV=1.027,R=1.25664 
INITIAL 
DERIVATIVE 
VT=R-C 
CDD=VT*KC 
CDOT=REALPL(0.0,1.0,CDD) 
C=INTGRL(0.0,CDOT) 
CONTRL FINTIM=0.01, DELT=0,.000001, DELS=0. 000002 
PRINT 0.0001, VT,CDOT,C 
SAMPLE 
CALL DRWG(1,1TINE,C) 
CALL DRWG(2,1,C,CDOT) 
TERMINAL 
CALL ENDRW(NPLOT) 
END 
STOP 


ip pe 





ITLE SSCOND ORDER DIGITAL POSITIONING SERVO 
STEP RESPONSE AND OUTPUT OF THE UP~DOWN COUNTER 


A IS A FACTOR USED TO ADJUST THE SIZE OF THE CURVE 
MULTIPLE RUNS 


INTHG RKSFX 
INTGER NPLOT, NCUR 
CONST NPLOT=1, NCUR=5,0P=1,RN=1,KV=1,0E7,A=20.0 
INITIAL 
KC=OP/6, 2831853 
DERIVATIVE 
BIT=C *KC 
IGBIT=AINT( BIT) 
NT=RN-IGBIT 
VT=NT/KC 
CDD=VT*KV 
CDOT=R=ZALPL(0.0,1.0,CDD) 
C=INTGRL(0,.0,CDOT) 
NO=NT*A 
CONTRL FINTIMH0.01, DELT=0.000001 , DELS=0.00002 
PRINT 0,0001,VT,CDOT,C, NT 
SAMPLE 
CALL DRWG(NPLOT,1,TIME, NO) 
CALL DRWG(NPLOT, 2, TIMZ,C) 
TERMINAL 
IF(NCUR. 52.1) GO TO 6 
NPLOT=NPLOT +1 
NCUR=NCUR-~1 
IF(NCUR. £Q.4) GO TO 
IF(NCUR, E2.3) GO TO 
DaecuRe).2) GO TO 
IF(NCUR,EQ2.1) GO TO 
1 OP=10,0 
RN=2.0 
A=1.5 
GO TO 5 
2 OP= 45.0 
Ri=9),0 
A=0,.3 
GO TO 5 
3 OP=90.0 
RN=18.0 
A-0.15 
GO TO 5 
4) OP=180.0 
RN=36, 0 
A=0.075 
5 CALL RERUN 
RETURN 
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»T 
EN 


nd 


STOP 


6 Cae 


NDRW(NPLOT) 


eas 





TITLE THIRD ORDER DIGITAL SERVO 

z VELOCITY AND OUTPUT OF UP-DOWN COUNTER 
* VARIATION OF SLOPE OF RAMP FUNCTION 

x 
* 


MULTIPLE RUNS 
* 
meas SKSFX 
INTGER NPLOT,NCUR 
CONST NPLOT=1, NCUR=4,KV=2. 5,0P=10.0, Be5.0 
Waa TA be 
KC=0P/6,.2831853 
DERIVATIVE 
RAM=B*RAlP(0.0) 
RN=AINT( RAM) 
BIT=C*KC 
IGBIT=AINT(BIT) 
NT=RN~IGBIT 
VT=NT/KC 
CD=VT*KV 
CDD=REALPL(0.0,0.2,CD) 
CDOT=REALPL(0.0,1,.0,CDD) 
C=INTGRL(0.0,CDOT) 
NO=N’ 
CONTRL FINTI*«M=30.0,DELT=0.004, DELS=0.03 
PRINT 0,2,VT,CDOT,C, NT 
SAKPLE 
CALL DRWG(NPLOT,1, TIME, NO) 
CALL DRWG(NPLOT, 2,TIME, CDOT) 
TERMINAL 
IF(NCUR, 2.1) GO To 6 
NCURS=NCUR=1 
NPLOT="PLOT+1 
IF(NCUR.£Q.3) GO TO 1 
IF(NCUR,EQ,2) GO TO 2 
IF(NCUR.EQ.1) GO TO 3 
1 B=4,0 
GO TO 5 
2 B=3.0 
GO TO 5 
Bis, 
5 CALL RERUN 
RETURN 
6 CALL ENDRW(NPLOT) 
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ITLE. THIRD ORDER POSITIONING SERVO 
STEP RESPONSE 


NULTIPLE RUNS 


*% %8 HR KR HHQ 


ANALOG AND DIGITAL SYSTEMS 
* 
INTSG RKSPX 
INTGER NPLOT, NCUR 
CONST NPLOT=4, NCUR=1,C=0,0, D=1,0,KV=6.0,RN=2.0,0P=10.0,P1=0.2,... 
Z1=0,5,P2=12,5,22=1,25 
INITIAL 
KC=0P/6, 2831853 
TZ1=1,0/Z1 
TP1=1.0/P1 
TZ2=1.0/Z2 
TP2=1.0/P2 
DERIVATIVES 
BIT=KC*C*D 
IGBIT=AINT( BIT) 
NT=RN-IGBIT-C*C 
VT=NT*(D/KC+C) 
VIN=VT*KV 
CD=LEDLAG(0.0,T21,TP1, VIN) 
C@R=LEDIAG(0,0, 722, TP2,CD) 
CDD=REALPL(0.0,0.2,COD) 
CDOT=REALPL(0.0,1,0,CDD) 
C=INTGRL(0.0,CDOT) 
CONTRL FINTI=20,.0,D=LT=0.0025, DELS=0.015 
PRUNE Oj2,VE,eDoT,C, NT 
SAMPLE 
CALL DRWG(1,NCUR,TIME,C) 
TERMINAL 
IF(KPLOT.&Q.1) GO To 6 
NPLOT=NPLOT~1 
NCUR=NCUR+1 
TF(NPLOT. EQ.3) GO 
IF(NPLOT, 22.2) GO 
IF(NPLOT, EQ.1) GO 
2 RNEI,0 
oPp=45.0 
GO TO 5 
3 RN=36.0 
OP=-180.0 
GO TO 5 
4 RN=1, 25664 
C=1,0 
D=0.0 
5 CALL RERUN 
RETURI 
6 CALL ENDR2(NPLOT) 


LOZ 
LO3 
TO 4 


END 
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